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ABSTRACT 

The recent discovery of super-Earths (masses < 10 M®) has initiated a dis- 
cussion about conditions for habitable worlds. Among these is the mode of con- 
vection, which influences a planet's thermal evolution and surface conditions. On 
Earth, plate tectonics has been proposed as a necessary condition for life. Here 
we show, that super-Earths will also have plate tectonics. We demonstrate that 
as planetary mass increases, the shear stress available to overcome resistance 
to plate motion increases while the plate thickness decreases, thereby enhanc- 
ing plate weakness. These effects contribute favorably to the subduction of the 
lithosphere, an essential component of plate tectonics. Moreover, uncertainties 
in achieving plate tectonics in the one earth-mass regime disappear as mass in- 
creases: super-Earths, even if dry, will exhibit plate tectonic behaviour. 

Subject headings: planetary systems — planets and satellites — Earth 
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Introduction 



Until recently, Earth was the largest terrestrial object known to exist. However, five 
super-Earth planets (a class defined as having a mass between 1-10 M^- earth-masses) have 



been detected in the las t few years (IRivera et al. II2005I : I Beaulieau et al. II2006I : iLovis et al 



20061 : lUdry et al. 1120071 ). The five planets have masses in the 5-10 range, but we do not 
have information on their sizes and cannot be sure if these are really rocky terrestrial planets. 
However, their discovery provides some evidence that super-Earths might be common and it 
is only a matter of chance that our Solar System has none. Some of these planets might be in 
the 'habitable zone', where the radiation from the star allows for the presence of liquid water, 
but only their thermal and chemical evolution will determine if they are, in fact, habitable. 
In turn, their thermal evolution and surface conditions depend on and affect their tectonic 
regime. Currently, Earth is the only planet where plate tectonics is active. Furthermore, 
this mode of convection has dominated our planet's geological history, is associated with 
geochemical cycles a nd thus, has been proposed as a required mechanism for life on Earth 
( iWalker et al. Ill98ll ). Here we address whether or not super-Earths are likely to have plate 
tectonics or be in a stagnant lid convection like Mercury and Mars. 



2. Analysis 

Plate tectonics is a complicated process that primarily requires lithospheric failure, de- 
formation and subduction. For this, convective stresses of the system have to be large enough 
to overcome lithospheric resistance to deformation and the plate has to reach negative buoy- 
ancy to drive subduction. The convective stresses are a function of the underlying fiow field 
and the viscosity of the fiuid, whereas the plate's strength depends on the elastic thickness 
and mechanical properties. Our contribution is based on identifying how the lithospheric 
or plate thickness (5) and convective stresses (r, a) depend on planet ary mass (M). ^y"e use 



detai led models of the internal structure of massive terrestrial planets (jValencia et al. 112006 



20071 ) and parameterized convection models to show that the condition for deformation and 



subsequent subduction is easily met by planets more massive than Earth. A second condi- 



( Sleep 


1992: 


Davies 


1993) 



requirement is with recent elaborated petrological models (jHynesll2005l : lAfonso et al. 1120071 ) 
suggesting that on Earth, the plate's mean density is never larger than the underlying man- 
tle's. Nontheless, we discuss the conditions for the development of negative buoyancy on 
super-Earths. 
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2.1. Parametric Convection 



The thickness of the lithosphere or boundary layer (6) depends on the Rayleigh number 
(Ra) - a parameter governing convection. A variety of mod els of boundary layer convection 
dO'Connell and Hager lEoSol : IXurcotte and Schubert Il2002h lead to 5 ~ D {Ra/Rac)\ with 
s = — 1/4, and Rqc the critical Rayleigh number above which a fluid starts convecting. This 
results in 5 being independent of the depth D of the convecting layer. 

The exponential relationship between 5 and Ra has been extensively addressed as dif- 
ferent flow details, ge ometries, turbulence , and other complicating effects are taken into 
account. In particular, IConrad and Hager I (1l999l ) suggested that s ~ where viscous dissi- 
pation within the bending subducting slab is lar ger than the dissipation wi thin the shearing 
mantle, which is the case for thick strong plates (jConrad and Hager Ill999al ). However, since 
we show below that the plates for super-Earths are thin, the case where s ~ can be ruled 
out. Other values of s will not change the qualitative results obtained here. 

The deviatoric horizontal normal stress (cr) responsible for causing failure on the plate 
is (to first order) balanced by the shear stress (r) applied over the base of the plate. Thus, 
by a simple force balance cr = r where L is the plate length, calculated from the time 
it takes for the plate to cool conductively and grow to a thickness 5, L ~ 5^ m/k, where 
u is the convective velocity. The shear stress depends on the viscosity rj and the velocity 
field of the fluid. In the most simple case this dependence will be linear so that the stress 
under the plate is t ^ rj u/D. The faster the flow, the larger the stress. In general these 
are competing effects: viscosity decreases while velo cities increase with higher tem peratures 
(T). The relationship for velocity isu ^ k/D Ra^^'^ ( Turcotte and Schubert II2002I ). A super- 
Earth can be expected to have larger velocities due to larger Rayleigh numbers. Thus, we 
need to investigate the extent to which the viscosity effect reduces the stress to predict 
whether or not plate failure will occur. We do so by considering two cases for viscosity: an 
isoviscous case (?7 = tjq = lO^^Pa-s - E arth's nominal value), and a T dependent viscosity 
case {t] = ?7o(T/To 



-30 



Davies I (Il980h l 



2.2. Structure and Convection 



W e calculate the interna l structure of terrestrial planets with masses in the range of 1 — 
IOM0 (jValencia et al. II2006I ) and determine the mantle thickness, density and gravity of each 
planet. We find that a power law relationship can adequately express the dependence of these 
parameters on M. For planetary mantles, the Rayleigh number depends on the mantle density 
(p = pe(M/Me)°-'°), gravity {g = g^{M/M^f'^), mantle thickness {D = D^{M/M^ 
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and heat flux (q), which all depend on the mass of the planet. Ra also depends on material 
properties: thermal expansivity [a), thermal diffusivity (/t), thermal conductivity (k) and 
viscosity [r]). The latter is a strong function of T. Therefore the Rayleigh number Ra = 
apgD^q/ ukr], depends indirectly on M and T. It increases for super-Earths as the size of their 
mantle and interior heat increases. To account for the heat fiux as a function of M we scale 
radioactive heat production and indirectly consider secular cooling. It is straightforward to 
scale the radioactive production with M for similar bulk compositions. Secular cooling, on 
the other hand, can only be completely assessed by considering the full thermal evolution. 
More massive planets would have larger initial gravitational energy budgets to dissipate over 
time; we assume for simplicity that the heat to be lost to secular cooling is proportional to 
M. Convcctivc velocities also scale roughly with M {u ~ M^'^^- see Table 1), so that the 
proportional rate of cooling should be roughly independent of M. Even though the relative 
contributions of radioactive heating and secular cooling to the total heat flow will depend 
on the evolution of the planet, we will assume that heat flow scales proportionately with 
M, especially because we expect the radioactive production to be dominant. The internal 
structure model allows us to calculate Ra, 6 and T beneath the plate for planets between 
1-10 M®. Convergence of these parameters determines the structure and final radius of each 
super-Earth planet. 



3. Results 

Figure 1 shows the dependence of plate thickness (blue) and horizontal normal stress 
(green) as a function of M for the T-dependent- viscosity case and Earth-hke radioactive heat 
sources (sohd). An additional case with reduced heat sources (dashed) is discussed below. 
The trends between these two scenarios and scaling exponents are similar, although their 
IM0 values are different. 

As expected, the plate thickness decreases as the M of the planet increases. This is 
because a more vigorous convective interior can transport heat more efficiently to the surface 
and sustain a higher surface heat flux. In addition, shear stress underneath the plate increases 
proportionately with M (Fig. 2), and thus, the dcviatoric horizontal normal stress also 
increases. This means that the viscosity reduction effect on stress is very small compared 
to the velocity effect. This is the case because the drop in T within the plate and hence 
local viscosity arc nearly independent of M as indicated by our calculations for the more 
realistic T-dependent- viscosity case (Fig. 2). For similar surface temperature, the increase 
in temperature at the base of the lithosphere between Earth and a super-Earth differs at 
most by only 5%. We think the reason for this surprising result is a very efficient negative 
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feedback induced by viscosity. If T increases, the viscosity decreases and heat flux increases, 
coohng the planet and reducing T. We ignore pressure (P) effects on viscosity because we 
expect them to be small under thin plates and because we find that P at the base of the 
plate is roughly the same for all super-Earths. 

Figure 2 shows the dependence of all the relevant convective parameters with plane- 
tary mass normalized to the Earth's values for the T-dependent case. The results can be 
adequately fitted to a power law such as 6 = 5(q{M/M^Y where the exponent [3 shows the 
dependence on M. In a vigorously convective interior (high Ra) the velocities are expected 
to be large and this is evident in Fig. 2. The relative size of the convection cells measured 
as the length of the plate with respect to its radius [L/R) increases slightly, from a 0.29 for 
a IMq planet to a 0.30 ratio for a lOM^ planet. The maximum plate age, related to the 
convective overturn, in a super-Earth is much less (see Table 1 for (3 = —0.91). 



3.1. Heat Flow 

In addition to a bulk silicate Earth composition for radiactive sources, we also consider 
planets with no "^^K content. As a result o f potassium's volatile nature. Earth is depleted in 



this element ( iMcDonough and Sun Ill995l ). It is therefore possible for other planets to also 
have non-chondritic K concentrations. U/0 and Th/0 ratios can be expected to be fairly 
constant for planets forming around stars with a relatively recent supernova event (within 
10 Gy). K contributes less than 1/3 of total heat from radioactive sources at the present 
time on Earth. Within this context, we consider a pesimistic scenario for a super-Earth with 
only 1/3 of its chondritic heat flow (Q) and no secular cooling contribution. 

This end- member case (Q/3) shows comparatively thicker plates (Fig. 1), with 88 km 
versus 43 km in thickness for a IM® planet. Additionally, the driving force is in the order of 
1-6 MPa which is lower than the Earth's value of ~10 MPa. This is an extreme case where 
subduction would be difficult but still possible for the most massive super-Earths if other 
factors, such as the presence of hydrated minerals is invoked to reduce the strength of the 
plate. For each planetary mass, there is a threshold in heat flow below which the convective 
driving force would not be sufficient to maintain subduction and plate tectonics would cease. 
Nontheless, most of the super-Earths that will be detected in the next few years will likely 
fall in a regime above the Q/3 scenario. 
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3.2. Lithospheric Buoyancy 



On Earth, the oceanic hthosphere is comprised of 7 km of basahic crust (chemically 
buoyant), overlaying lithospheric mantle. We perform a simple buoyancy calculation at sub- 
duction ages by considering the thermal contraction from the ridge. We determine maximum 
crustal thicknesses that still allow the lithospheric mean density to be larger than the underly- 
ing mantle's. We find this value to be 13% the plate thickness at subduction zones (twice the 
average thickness - 6). The crustal thickness depends on the extent of melting under ridges 
which depends on the P and T conditions. As shown above, T at ridges and under plates 
varies little with M, so that with larger gravity values, melting occurs at shallower depths in 
super-Earths. A simple calculation shows that the crust thins with increasing planetary M, 
and while it becomes a larger component of the also thinning boundary layers, its thickness 
stays below the threshold allowing negative buoyancy. This effect stands in contrast to the 



young Earth that had a higher potential T but same P profile (INisbet and Fowler I Il983 



Davies 



those of 



9931). To pre cisely assess the den s ity of the plate detailed petrological models hke 



Hyned (120051 ) and lAfonso et al 



(120071 ) are needed. It is important to point out 
that their models do not predict negative buoyancy on Earth when the plate's density is 
compared to sublithospheric density, arg uing that perhaps thi s condition is not critical for 
the initiation of subduction. In addition, iBecker et al. I (119991 ) showed that in compressive 
regions, the cold plate can thicken to a point that negative buoyancy is inevitable and sub- 
duction follows. If negative buoyancy is indeed a necessary condition, our simple calculations 
show that super-Earths can satisfy it. 



4. Summary and Implications 

In summary, convection is more vigourous in massive terrestrial planets, making their 
lithospheres thinner and therefore reducing lithospheric strength. Furthermore, they achieve 
larger stresses owing primarily to larger velocities and therefore can more easily overcome 
the lithospheric resistance to deformation. Plates may reach negative buoyancy on super- 
Earths despite their relative younger ages. This scenario is suitable for the failure of the 
plate and subsequent subduction, which is a necessary step for plate tectonics. Given that 
Earth's convective state leads to plate tectonics, the more favorable conditions experienced 
by super-Earths will inevitably lead to plate tectonics. Furthermore, planets of similar mass 
should have the same potential to exhibit plate tectonics. Conversely, this physics can help 
explain why small planets like Mars, Mercury and the Moon do not exhibit plate tectonics. 
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4.1. Role of Water 



Venus is o nly slightly smaller than Earth and d oes not exhibit plate tectonics, although 
some authors (iTurcotte Ill993l : iJellenik et al. II2005I ) have suggested it may have in the past. 
This observation indicates that the ~1 earth-mass case falls within a zone of transition 
between 'hard' stagnant lid and mobile plate regimes. In this case, characteristics other 
than M may be important to the dynamics of the lithosphere. For example, the high surface 
temperature of Venus might lead to a weak, highly deformable boundary layer that would 
not support the coherent plate-like behaviour that characterizes oceanic plates on Earth. 
Moreover, plate strength is relatively large compared to the mantle driving force in the one 
earth-mass case; yield stresses are on th e order of 1-5 GPa for olivine (the representative 
upper mantle mi neral) (IChen et al. II19981). wh ereas our calculations, in agreement with more 
detailed models (IBecker and O'Connell Il200ll ). suggest an underlying driving force of only 10 
MPa. Since slip can occur on pre-existing faults at stress values of a few MPa, the existence 
of plate tectonics in the one earth-mass regime may thus depend crucially on the conditions 
required to initiate subduction. The presence of water is one possible mechanism to reduce 
the yield strength of a plate and friction on faults. Experiments show that water reduces 
the yield strength of olivine by 62% when raisin g the temperature from 25 to 400°C at 10 
GPa, compared to a drop of 39% in dry olivine (jChen et al. Ill998l ). Hence, the hydration 



level of Venus' mantle, which is 1-2 orders of magnitude lower than on Earth (jZolotov et al. 



19971 ). may make it very difficult for convective forces within this planet to overcome plate 
resistance. 



For larger planets, super-Earths, these issues become less relevant. A wet super- Earth 
will clearly have enough driving force to sustain subduction. But, more importantly, the 
consequences for initiating subduction associated with the hydration of a one earth-mass 
planet (i.e., a reduction of the yield strength by half) would be similar to a doubling of the 
mass of the planet (Fig. 2). That is to say, both scenarios would be as likely to initiate and 
maintain subduction. 



4.2. Atmospheric Observables 



The difference between a Super-Earth with active plate tectonics and one with stagnant 
lid is in the access of upper mantle material and gasses to the atmosphere. The ffist case 
allows several global geochemical cycles to operate, like the CO2 and SO2 ones. For example, 
cases in our Solar System compr ise: Earth with a CO2 cycle and possibly early Mars with 
a SO2 cycle ( Halevy et al. 2007h. Ear t h has had stable modest levels of atmospheric CO2 
(between 160-7000 ppm - iRoyer et al. I (I2OOII )) in the last 0.5 Gy whereas Venus levels stand 
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today at 96%. A planet with plate tectonism and a carbonate rock reservoir has a n efficient 



built -in cycle that stabilizes climate at temperatures within the liquid water regime (IKasting 



19961 ). A super-Earth that has plate tectonics and weathering capabilities can be expected 
to have CO2 atmospheric concentrations that would yield temperatures around liquid water. 
Therefore, evidence against the presence of plate tectonics on an exoplanet would be the 
detection of high values of CO2 for the age of the star, type of star and orbital distance. An 
SO2 based atmosphere is also possible and the same reasoning would apply, since the sulfur 
cycle operates analogously to the carbon cycle. But obviously, more theoretical research is 
necessary to model the details and predict the right observable signatures. 

In conclusion, we show here that as mass increases, the process of subduction, and hence 
plate tectonics, becomes easier. Therefore, massive super-Earths will very likely exhibit plate 
tectonics. In the future with TPF by NASA and Darwin by ESA it might be possible to use 
spectroscopy to identify atmospheric signatures suggesting plate tectonism on these objects. 
This class of planets offers the possibility of finding Earth analogs and, in particular, make 
attractive targets in the search for habitable planets. 
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Fig. 1. — Plate tliickness and stress for super-Eartlis. Tlie dependence of plate thickness 
(blue and left axes) and horizontal normal stress (green and right axes) on mass are shown 
for a BSE composition (solid) and a reduced heat content case (dashed). Both cases have 
a T-dependent viscosity treatment. As mass increases, the stress available to overcome 
plate resistance increases, while plate thickness decreases, enhancing pl ate weakness. Plate 



2000 



thickness for Mars and Venu s is shown with their uncertainty (from (IZuber et al 

Nimmo and McKenzie I [igosl )). The Earth's plate thickness is a glob al average. Planet s 
with masses larger than lOM^ can start to retain gas during formation (jida and Lin 112004 ). 
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Fig. 2. — Convective parameters scaling with planetary mass. Plate thickness (solid line), 
drop in temperature within the plate (dash-dotted line), shear stress (medium dashed line), 
plate length (long dashed line), normal stress on plate boundaries (dotted line) and convective 
velocity (short dashed line) depend on the mass of the planet. This dependence can be fitted 
to a power law with corresponding exponents shown as numbers. The values are normalized 
for the Earth - see Table 1. 
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Table 1. Convective parameters' power law dependence on mass 



Parameter 


Earth 


77 = 77(r) 


Isoviscous 


plate thickness (5) 


43 km 


-0.45 


-0.29 


drop in temperature (AT) 


1277 °C 


0.02 


0.18 


plate's velocity (u) 


3 cm/y 


1.19 


0.87 


shear stress (r) 


0.3 MPa 


0.27 


0.58 


normal stress (a) 


10 MPa 


1.00 


1.16 


plate length (L) 


1800 km 


0.28 


0.29 


convective time {t) 


70 myr 


-0.91 


-0.58 



Note. — The second column displays Earth's convective pa- 
rameters. The only assumed value is the plate velocity. The 
third and fourth column show the value for the exponent f3 in 
the power law relationship - ie. 6 = 5q){M / M^Y - for a BSE 
composition with a T-dependent viscosity and isoviscous cases 
respectively. 



